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Abstract 

The matter effects for solar neutrino oscillations are studied in a general scheme 
with an arbitrary number of sterile neutrinos, without any constraint on the mixing, 
assuming only a realistic hierarchy of neutrino squared-mass differences in which 
the smallest squared-mass difference is effective in solar neutrino oscillations. The 
validity of the analytic results are illustrated with a numerical solution of the evo- 
lution equation in three examples of the possible mixing matrix in the simplest case 
of four-neutrino mixing. 
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1 Introduction 



The possible existence of sterile neutrinos [1] is one of the most interesting open ques- 
tions in neutrino physics (see Refs. [2-8]), which may give important information on 
the physics beyond the Standard Model (see Refs. [9, 10]). Sterile neutrinos are neutral 
fermions which do not interact weakly, as the three ordinary active neutrinos z/g, t'^t, 
do. Moreover, sterile neutrinos can mix with active neutrinos, leading to active-sterile 
neutrino oscillations which are observable either through the disappearance of active neu- 
trinos or through flavor transitions generated by a squared-mass difference Am^ which is 
larger than the solar (SOL) and atmospheric (ATM) mass differences (see Refs. [2-8]) 

AmloL-8xlO-^eV^ (1) 
Am^TM - 2.5 X 10-3 eV^ . (2) 

In fact, the oscillation data show that z/g, t'^, Vt are mainly mixed with three light 
neutrinos z/i, z/2, z^s with masses such that 

^"^SOL = ^°^\\ . (3) 

^"^ATM = I ^"^31 1 - I ^"^32 1 > (4) 

with ^mlj = ml—m'j. The order of magnitude of the absolute values of the three neutrino 
masses mi, 1712, is bound to be smaller than about 1 eV by /3-decay, neutrinoless 
double-/? decay and cosmological data (see Refs. [6-8,11]). 

The possible existence of sterile neutrinos has been discussed by many authors in 
connection with the indication of — > Pg oscillations with Am^ ^0.1 eV^ found in the 
LSND experiment [12]. Since this squared-mass difference is much larger than Am^rj^j^j, 
the LSND oscillations can only be explained with mixing of more than three neutrinos, 
leading to the unavoidable introduction of additional sterile neutrinos. 

The LSND signal has not been seen in other experiments and is currently disfavored 
by the negative results of the KARMEN [13] and MiniBooNE [14] experiments [15-17]. 
On the other hand, indications in favor of active-sterile neutrino oscillations [18-21] come 
from the anomalous ratio of measured and predicted ''^Ge observed in the Gallium ra- 
dioactive source experiments GALLEX [22, 23] and SAGE [24-27] and from the Mini- 
BooNE [14] low-energy anomaly. Moreover, sterile neutrinos may constitute the cosmo- 
logical dark matter and active-sterile oscillations may explain the pulsar velocities [28-35] . 
Since the standard cosmological scenario constrains the masses, mixing, and number of 
neutrinos [?,?, 38-42], sterile neutrinos may be heralds of alternative cosmological mod- 
els [43-45]. 

Bearing in mind that the existence of sterile neutrinos which are not related to the 
LSND anomaly remains possible, the quest for sterile neutrinos is one of the best tools 
for the exploration of physics beyond the Standard Model. A convenient way to re- 
veal the existence of sterile neutrinos is the search for disappearance of active neutrinos 
due to active-sterile oscillations [46, 47] . The explanation of the Gallium radioactive 
source experiments anomaly and the MiniBooNE low-energy anomaly through active- 
sterile neutrino oscillations requires the disappearance of electron neutrinos into sterile 
states [18-21]. If these transitions exist, they contribute also to the disappearance of 
solar electron neutrinos. This effect, which has not been studied in detail so far, is the 
object of study of this paper. 
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The matter effects for solar neutrino oscillations in a four-neutrino scheme was studied 
in Ref. [48], where it was assumed that the electron neutrino has non- negligible mixing 
with only two massive neutrinos, z/i and z/2, which generate the solar squared-mass dif- 
ference AmgQL = Amli- In this approximation, electron neutrinos cannot oscillate into 
sterile states in the short distances required for the explanation of the Gallium radioactive 
source experiments anomaly and the MiniBooNE low-energy anomaly. Hence, a com- 
bined analysis of the Gallium radioactive source experiments anomaly, the MiniBooNE 
low-energy anomaly and solar neutrino data in terms of active-sterile oscillations [49] 
requires an extension of the study in Ref. [48] to the general case in which there is no 
constraint on the mixing of the electron neutrino. 

In this paper we study the matter effects for solar neutrino oscillations in a general 
scheme with an arbitrary number of sterile neutrinos, without any constraint on the 
mixing. Hence, we consider the general case of mixing of three active neutrino fields (z/g, 
u^, Vt-) and A^^ sterile neutrino fields: 

N 

l^aL = '^UakJ^kL {a = 6, fl, T, Si, . . . , SnJ , (5) 

k=l 

where N = 3 + Ng, Uk is the field of a massive neutrino with mass m^, and the subscript 
L denotes the left-handed chiral component. The numbering of the massive neutrinos is 
chosen in order to satisfy Eqs. ([3j) and (jlj) with the first three neutrino masses and we 
assume that the additional neutrino masses are much larger than mi, m2, ms, in order 
to generate short-baseline oscillations which could explain the Gallium radioactive source 
experiments anomaly and/or the MiniBooNE low-energy anomaly [18-21] and could be 
observed in future laboratory experiments [15, 46, 47, 50-55] and through astrophysical 
measurements [56-60]. Hence, we have the following hierarchy of squared-mass differ- 
encesQ 

Am2i < lAmgJ < Aml^ for k>3. (6) 

The N X N unitary mixing matrix U can be parameterized in terms of 3 + SNg 
mixing angles (see Ref. [7]). Three mixing angles mix the three active neutrinos among 
themselves and the other SNg mixing angles mix the three active neutrinos with the A^^ 
sterile neutrinos. There are also 1 + 2Ns Dirac phases which affect neutrino oscillations 
and 2 + Ng Majorana phases which do not affect neutrino oscillations [62-64]. However, 
in this paper we do not need an explicit parameterization of U. We derive the flavor 
transition probabilities of solar neutrinos in terms of the values of the elements of the 
mixing matrix. 

The plan of the paper is as follows. Section [2] is the main one, where we derive 
the flavor transition probabilities. In Section [3] we discuss the peculiarities of the flavor 
transition probabilities in the extreme non-adiabatic limit. In Section H] we illustrate 
the validity of our results in the simplest case of four-neutrino mixing, considering some 
examples in which we compare the values of the analytic approximation of the flavor 
transition probabilities derived in Section [2] with those obtained with a numerical solution 
of the flavor evolution equation. Conclusions are drawn in Section [51 

^ The different case of mixing of two active neutrinos and a sterile neutrino with two smaU squared- 
mass differences, both of which contribute to solar neutrino oscillations, has been studied in Ref. [61]. 
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2 Evolution of Neutrino Flavors 

Solar neutrinos are described by the state 

= X] i^c{x)\Va) : (7) 



o=e,/^,T,si,...,sjv. 



where x is the distance from the production point in the core of the sun and 

^a(0)=5<,e, (8) 

with the obvious normahzation = 1- 

From the mixing relation ([5]) between flavor and massive neutrino fields and the 
structure of the weak charged-current Lagrangian, 



N 



Ccc = -^ J2 5Z4L7''f/afe'^fcLW^; + H.c., (9) 

^ a=e,^,T k=l 

it follows that in vacuum the fiavor states {ua) which describe neutrinos created in 
charged-current weak interactions are unitary superpositions of massive states |z/^) which 
are the quanta of the corresponding massive neutrino fields (see Ref. [7]): 

N 

K)=^KkWk) {a = e,n,T,si,...,SNj. (10) 
fc=i 

The evolution of the fiavor transition amplitudes ipa{x) is given by the Mikheev- 
Smirnov-Wolfenstein (MSW) equatioiJl [65,66] (see Ref. [7]) 

4* = ^(C/A.=£/t+^)*, (11) 

where E is the neutrino energy and 

= {1pe,i^^„'^pT,1ps^,■■■,'ipsJ,y , (12) 
= diag(0, Amj-^, Am^^, Aml^, Am%^) , (13) 
A = diag(Acc + ^Nc, ^Nc, ^Nc, 0, . . .) , (14) 

with Amlj = ml — m| and 

Acc = 2EVcc , ^Nc = '2EV^c ■ (15) 
The charge-current and neutral-current matter potentials are given by 

^cc = V2GfN, ~ 7.63 X 10"^^ tt^^ eV , V^c = -^v^G'FiV„ , (16) 



We omit a diagonal term E + m\/2E which generates an irrelevant common phase. 
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where Gp is the Fermi constant, Ne is the electron number density, is the neutron 
number density, and A^a is Avogadro's number. It is convenient to use the standard 
definition of the electron fraction 

In a neutral medium we have 

iVe = -iVAn, Nn = -NA{1-Y,) . (18) 

g g 

where p is the density. From Eqs. ([T6|) and f[T8l) . we obtain 

^Nc = RncAcc , with i?Nc = ■ (19) 

For solar densities we have 

^cc ~ I^NcI ~ Am^^ < |Am^J for A; > 3 . (20) 

In order to decouple the fiavor transitions generated by Aqc ~ I^ncI ~ ^f^2i from those 
generated by the larger squared-mass differences, it is useful to work in the vacuum mass 
basis 

={^l...,^lf = U^<I>, (21) 
which satisfies the evolution equation 

^ Avi/V = ±.^M^ + UUU) vl/V , (22) 
dx zii/ 

In this basis, the solar neutrino state in Eq. ([7]) is given by 

N 

wi^)) = Y.'^nxM), (23) 
fc=i 

In the vacuum basis the inequality in Eq. fl2Ul) implies that the evolution of each of 
the amplitudes ipj, . . . , ipjf is decoupled from the others, since the coupling due to the 
matter effects is negligible. Therefore, we have 

i^^ki^) ^ ^kiO) exp (^-z , for > 3 , (24) 

with 

^.'^(0) = U:, , (25) 

from Eqs. §^ and 

On the other hand, the evolution of ipY and ip2 is coupled by the matter effect. In 
order to solve it, we consider the truncated 1-2 sector of Eq. (l22l) : 



'dx VEa U^iUl^Ac. Ami, + l^^a2|M J 
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V ) , (26) 



where the sums over the flavor index a cover only the active neutrinos, which have a 
matter potential in Eq. f|T^ . It is convenient to subtract the diagonal term 



'^"^ 



21 



(27) 



which generates an irrelevant common phase. The resulting effective evolution equation 
can be written as 



'dx W 



1 /-Am^^ + Acos2^ 
AE 



^21 

A sin 2^e-*^ 



A sin 2^e*'^ 

''21 



Am2i — v4cos2,^ 



V 



with 



where 



AcOs2e= 5^(|t/air-|f/a2r)^a = AccX, 

A sin 2^6^^ = 2 KiUa2A^ = AccYe''' , 

a 

X=|f/elp-|f/e2r + i?NC ( | t^al | ' " | f/a2 | 

Ns 

= |f/elP - |f/e2p - RnC ^ (|f/sal' " \Us.2\') , 



(28) 

(29) 
(30) 



(31) 



1=1 



Y 



U*iUe2 + RnC ^ U*iUa2 
U*iUe2 - -RnC ^ U*^lUs,2 



i=l 



(32) 



In these expressions we used the unitarity of the mixing matrix in order to express the 
neutral-current contribution either through the active neutrino elements of the mixing 
matrix or through the sterile neutrino elements. 

One can immediately notice that, apart from the complex phase e**^, the effective 
evolution equation fl28l) has the same structure as the effective evolution equation in the 
vacuum mass basis for Ue-i^ti or i^e-^r two-neutrino mixing with the two-neutrino mixing 
angle replaced by the effective mixing angle ^ and Aqc replaced by A, with 



Y 

tan 2^ = — 
X 



A = AccVX^ + Y^^ 



(33) 



In fact, one can easily check that in the limit of z/g-z/^ or i^e-'^r two-neutrino mixing, ^ 
coincides with the two- neutrino mixing angle and A becomes Aqq- Therefore, the effective 
mixing angle ^ plays the same role as the usual two-neutrino mixing angle in the solution 
of the flavor evolution equation. However, one must keep in mind that in general ^ is not 
a pure mixing angle given by a specific parameterization of the mixing matrix, because 
it depends on the neutral-current to charged- current ratio -Rnc- 
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For simplicity, we consider a medium with constant electron fraction Ye, which implies 
a constant neutral-current to charged- current ratio i?NC- In this case, the effective mixing 
angle ^ and the phase if are constant. Moreover, we can use the phase freedom in 
the parameterization of the mixing matrix in order to eliminate (p. In fact, the flavor 
transition probabilities obtained from the MSW evolution equation flTTl) are invariant 
under the phase transformations 

U^^ ^ e'^-U^^e''^" . (34) 
which gives, from the definition in Eq. fl30l) . 

A sin 2^e*'^ A sin 2^6'^'^+'^^-'^'^ . (35) 
Therefore, choosing ipi — ip2 = V^, the evolution equation (128|) is simplified to 

■A. /^^n = 1. (-^^li + ^cos2e Asin2e \ . X 

'fix V A sin 2^ Am2, -ylcos2ey' ' 

which has the same structure as the effective evolution equation in the vacuum mass basis 
for Ve-Vfj, or Ve-^T two-neutrino mixing. 

In practice. Ye and i?NC are not constant in the Sun, as shown in Fig. [TJ The electron 
fraction is almost constant in the radiative (0.25 ^ ^ 0.7) and convective (r > 0.7) 
zones where the Hydrogen/Helium mass fraction is practically equal to the primordial 
value. In the core (r < 0.25) the electron fraction increases towards the center because 
of the accumulated thermonuclear production of Helium during the 4.5 Gy of the life of 
the Sun. However, it is plausible that the effects on flavor transitions of the variation 
of i?NC in the core is negligible if the transitions occur mainly in a resonance which is 
located in the radiative or in the convective zone. Hence we solve the evolution equation 
of neutrino flavors in the approximation of a constant -Rnc- In Section H] we will show 
with numerical examples that the effects of the variation of i?NC in the Sun are indeed 
negligible. 

The effective Hamiltonian in the evolution equation fl36l) is diagonalized by the rota- 
tion 

_ ( cosu; sinci;\ fifj^'^ 
J \ — sina; cosu J \V'2^ 



(37) 



with 



J-^ Dili z-c, , 

tan2u; = "2 — ~o7 • 



A sin 2^ 
Amli — A cos 2,^ 

The angle u is the mixing angle between the vacuum mass basis and the effective mass 
basis in matter. Its expression is similar to the one in t'g-^'/x ^e-^^r two-neutrino mixing. 
It is therefore useful to note that the evolution equation corresponds to the evolution 
of the effective interaction amplitudes 



cos^ sin A 



(39) 



for which the matter potential is diagonal: 



^ (-Arnl^ COS 2i + A Am^, sin 2^ \ ( i'X^ ...x 
dxW2 AE\ Ami, sin 2e Am^, cos 2^ - A i ' ' ^ ^ 
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Figure 1: Hydrogen mass fraction Xh, electron number fraction Ye = (1 + Xh)/2, 
and neutral-current to charged-current ratio i?NC [Eq. ffT^ ] in the BP04 Standard Solar 
Model [67]. -Rsun ^ 6.7 x 10^° cm is the radius of the Sun. 



This basis of effective interaction amplitudes is physically important for the study of the 
flavor evolution, because it is the basis in which the mixing generates transitions through 
the off-diagonal term in the evolution equation. Moreover, in this basis one can see that 
there is a resonance when the diagonal terms are equal, i.e. for A = A-^, with 



= Amgi cos 2,^ . 



(41) 



The relations between the effective mass and interaction bases in matter and the flavor 
basis are discussed in Appendix Rl 

The connection between the effective interaction amplitudes ipl and ipl and the am- 
plitudes ipf^ and of the effective energy eigenstates in matter is similar to that in 
two-neutrino mixing: 

^ f cos^M sin^M\ fi'Y^ 
where is the effective mixing angle in matter given by 



(42) 



From Eq. (!38|) we obtain 



tan 2^: 



M 



tan 2^ 



1 - (A/Aml^ cos 20 



(43) 



(44) 
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The effective mixing is maximal (^m = ^/4) in the resonance A = A^. Notice that 
Eq. (jH]) and the resonance condition fj4T|) are similar to the corresponding equations 
in z/e-t'^ or i^e-^r two- neutrino mixing (see Ref. [7]), with the two- neutrino mixing angle 
replaced by the effective mixing angle ^ and Aqc replaced by A. 

The interesting limiting cases for the effective mixing angle in matter are: 

^M^e, ^^0, for A<Ar, (45) 

eM^|, for A>Ar. (46) 

The evolution of the amplitudes ip^ and ip^ of the effective energy eigenstates in 
matter is given by 

dx [i'fj ~ AE [UEdcu/dx Ami, J [i^^ J ' ^ ' 

with the effective squared-mass difference 



Amlj = y (Am^i cos 2^ - A) + (Am^^ sin 2^ • (4. 
The non-adiabatic off-diagonal transitions depend on 

du Aml^A sin 2^ d In Ne 



- (49) 
dx 2(Am2^)' dx ' 

It is convenient to define the adiabaticity parameter 

^ Ami, ^ .5Q^ 
' 4:E\duj/dx\ 2EAm2^Asin2^|dlnA/'e/dx| ■ ^ ' 

The evolution is adiabatic if 7 ^ 1 all along the neutrino path. 

The solution of the evolution equation fHTl) leads to the averaged probabilities 



K(x)p = iv^f (o)r (1 - P12) + li^mr P12 , (51) 
Wi^ = iv^f (o)p P12 + li^fm' (1 - P12) , (52) 

where P12 is the probability of non-adiabatic off-diagonal transitions in Eq. (1471) . The 
averaging process washes out the interference between ip^{x) and ip^{x), which is not 
measurable because of the energy resolution of the detector and the uncertainty of the 
production region. The initial values (0) and are given by Eqs. (I25p and (13 7p 

with the effective mixing angle in the production region: 

(0) = coscuV^(O) -sincjVy(O) = cos f/*i - sin f/^g , (53) 
ipf{0) = sinwV^(O) + coscjVy(O) = sin uj^U*^ + coscu°f/*2 . (54) 

We can use the phase freedom in Eq. (IMI) in order to have real Uei and 11^2, keeping the 
reality of 2 J2a f^ai^a2^a in Eq. (150]) . since there are three arbitrary phases available: ipe, 
ipi, and (p2- With this choice, Uei and Ue2 can be written as 

N 

Uei = COS^e COSXe , ^e2 = siu COS Xe , with siu^ Xe = ^ |f^efcP • (55) 

k=3 



Then, ip^{0) and ip^{0) are given by 

(0) = COS?9°COSXe, ^^(0) = sint9°cosXe, 
with the effective mixing angle in the production region 

From Eq. (I97|) . the electron neutrino state is given by 



(56) 



(57) 



N 



= [cos(^9e + uj)\vf) + sin(?9e + u:)\vf)] cosXe + ^*c.kWk) 

k=3 

N 

= [C0s(^e - owl) + M^e - owl)] COSXe + K.^^) . 



{51 



fc=3 



Hence, an electron neutrino contains both |z/}) and \i>l) if ^ 7^ "^g. In general, this is 
the case, since the effective mixing angle ^ in the effective evolution equation (125]) of the 
truncated 1-2 sector is different from the electron neutrino mixing angle -^e because of the 
neutral-current contribution to the matter potential (see the discussion after Eq. fITOl) ). In 
fact, as one can see from Eqs. fl3Tl) - fl33l) . the effective mixing angle ^ does not depend only 
on the mixing of u^, but also on the mixings of z/^ and u^- which affect the neutral-current 
contribution to the matter potential. 

Since the matter effects in the detector are negligible, the averaged flavor transition 
probabilities are given by 




N 



k=l k=3 
2 N 

E WWW' + E W\^ki^)\ 



k=l 
2 



k=3 

N 



(59) 



k=l k=3 

where we have used Eqs. flM|l and fl25l) . Writing jt/gip and |f//32p as 

N 

|?7/3i|^ = cos^t9/3 cos^X/3, |?7/32|^ = sin^^9/3 cos^X/3 , with sin^ X/3 = E '^^^'^ ' ^^'^^ 



k=3 



and using Eqs. (15Ti) . and (IS^ . we obtain 



+ 



P12 1 cos2i9flCos2?9° 



N 



COS^ Xl3 COS^ Xe + E I ^/^fc n t^efc I 



(61) 



k=3 
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It is clear that in the hmit of two- neutrino mixing, in which cos^Xe = cos^X/j = 1, the 
transition probabihty reduces to the well-known Parke formula [68] (see Ref. [7]). 

For neutrinos produced above the resonance, where A ^ Ar, from Eqs. fH6|) and fl57|) 
we have 

+ ^ cos2^9° ~ -cos2(^9e-0- (62) 



In this case. 



(\ \ 1 ^ 

- ( --P12J COs2^9;3COs2(t9e-0 COS^X/SCOS^Xe + JZltZ/jfcnf/.fcp. (63) 
^ ' -I fc=3 

The crossing probability Pyi can be inferred by using the similarity of the evolution 
equation (H7I) with that in two- neutrino mixing [69-72] (see Ref. [7]): 



exp (-f 7rF) - exp [-\ 
1-exp (-f7R-^ 



F 



P12 = 7 ^ X ^(^0 - Ar) , (64) 



where 7r is the adiabaticity parameter at the resonance (HTl) . which corresponds to the 
minimum of Am^^: 

(Ami^)j^ = Am^isin2e. (65) 
Then, from Eq. (150|) . 7r is given by 

^ Am^ sin^ 2^ 

2Ecos2e|dlnArjdx|R' ^ ^ 

The value of the parameter F depends on the electron density profile. For an exponential 
density profile, which is a good approximation for the solar neutrinos, [69-77] 

F = 1 - tan^ i . (67) 

In Eq. fl64l) there is a 6'-factor which reduces P12 to zero when the neutrino is created at a 
density which is smaller than the resonance density {Aq is the value of A at production). 

Let us finally notice that for solar neutrinos it is sufficient to calculate the v^. sur- 
vival probability P^^ and the sum of the transition probabilities into sterile neutrinos 
ue^fa-: from which Y2f3=^,T P'^e^'^/} trivially obtained using the conservation of 
probability. The separate transition probabilities into z/^ and Ur are not needed because i/^ 
and u-r have the same interactions in the detector. In fact, and Ur can be distinguished 
only through charged-current interactions with production of the associated charged lep- 
ton, which is forbidden by the low energy of solar neutrinos {E < IbMeV). Hence, z/^ 
and Ur are detected only through neutral-current interactions which are flavor-blind. 



3 Extreme Non-Adiabatic Limit 

It is interesting to note that in the extreme non-adiabatic limit 7r -C 1 the transition 
probability Pu^^y^^ is different from the averaged transition probability in vacuum (VAC) 



N 



k=l 



11 



1 COS 2{}g COS 2'(9e 

2 2 ^ 



N 



COS^ X/3 COS^ Xe + ^ l^^/^fcHf/efcl 



fe=3 



In fact, since in the extreme non-adiabatic limit 



(7R«1) 
12 



COS ^ - 



from Eq. (!63|) we obtain 



-p(7R«l) _ 



^ + i COS cos 2(?9e - COS 2^ 



TV 



COS^ X/3 COS^ Xe + ^\ U(3k\^\Uek\ 
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(69) 



(70) 



k=3 



This expression coincides with that in Eq. fl68p only if ^ = -i^g- From Eqs. fl31l) - fl33p 
one can see that ^ = "i^g if there is no neutral-current contribution in Eqs. fl31l) and fl32p 
or if the neutral-current contributions cancel in the ratio Y/X. The first possibility is 
obviously realized in a neutron-free medium (-Rnc = 0). It is also realized if Usi = 

Us2 = for S = Si,...,Sn,, which implies that Z]a=e,a,r l^^lP = I]Q=e,Ai,r l^"2p = 1 



and Yla=e ^iT^aiUa2 = 0. This is obviously the case in three-neutrino mixing (no sterile 
neutrinos). The second possibility is realized in four-neutrino mixing with [/es = f^e4 = 
[48]: in that case 



Uel = COS-de 



e2 



sin 'de 



and 



f/^l = - sin ■(9e cos Xs , Us2 = COS'&eCOSXs ■ 

Using the unitarity of the mixing matrix, we have 

J2 - If^«2|') = \Us2\' - \Usl\' = COs2^eCOs2xs 

2 Yl KiUa2 = ~U:,Us2 = sm2'deCos\s, 



(71) 
(72) 

(73) 
(74) 



leading to 



X = cos 2'de (l + -Rnc cos^ Xs) 

Y = Sm2l^e (1 + i?NC cos^ Xs) 



(75) 
(76) 



and tan 2^ = X/Y = t&n2i9e. 

In the case of two-neutrino mixing the averaged flavor transition probability in the 
extreme non-adiabatic case coincides with the averaged flavor transition probability in 
vacuum because the flavor states are conserved in the resonance. A z/g created above the 
resonance is practically equal to a z/^, which travels undisturbed to the resonance, where 
it is still a (it only develops a harmless phase). If the crossing of the resonance is 
extremely non-adiabatic, the z/g emerges unchanged from the resonance (apart from an 
irrelevant phase) and the flavor transition measured on Earth is only due to the vacuum 
oscillations from the Sun to the Earth. On the other hand, in the case of mixing with 
sterile neutrinos, we have seen in Eq. (158|) that, if ^ 7^ "(^e, a ^'e contains both z/{ and 
i>l, which are conserved if the crossing of the resonance is extremely non-adiabatic. The 
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different phases developed by z/{ and before and during resonance crossing make the 
averaged flavor transition probability different from that in vacuum. We can show that 
explicitly following the evolution of the solar neutrino state \i>{x)). If the neutrino is 
created clS db i^g above the resonance, the initial state is, from Eqs. fHBl) . fISS]) . and (1971) 



N 



Wm = We) = [- Sin(^e - OK) + COs(^e - OW^)] COS Xe + t/^.I^J) • (77) 



k=3 



Since from the production to the resonance the effective massive states develop different 
phases, just before the resonance the solar neutrino state is 



\u{x^)) = [-sin(?9e - Oe''^'^"^Vf ) + cos(?9e - ^e^'^'^'^^W^) 



COSXe 



N 



k=3 



Since above the resonance |z/{) ~ |z/^) and {uD ^ — |z/^), we have 

= [sin(79e - Oe"^'^"^V2) + cosine ' ^e^'^'^'^Vl) 

N 

fc=3 



COSXe 



(79) 



If the crossing of the resonance is extremely non-adiabatic, there are no transitions among 
the effective interaction states |z/{) and l^l), which develop only different phases. Hence, 



just after the resonance we have 

;+)) = [sin(t9e - Oe'^'^''^V2) + cos(?9e - Oe' 

N 



UiXr 



i<^2(x+)| I 



COSXe 



^0) 



A:=3 



Since after the resonance the neutrino propagates practically in vacuum, it is convenient 
to express its state in terms of the propagating massive states: 

\u{x > xr)) = [- sin(^9e - sin^e^^i^^') + cos(^9e - cos^e^'^^^^)] cosxe^i^) 
+ [sin(^9e - cos^e^'^i^^) + cos(^9e - sin^e^'^^^^)] cosXe^^) 



AT 



^1) 



A:=3 



Since 0i(x) ^ 02 (a;), if ^ 7^ "i^e the absolute values of the amplitudes of lu^) and jz/g ) are 
different from those of 

N 

I u,) = [cos I ^^i"") + sin A I ^2)] COS Xe + 5^ U*^, I U^) . (82) 



fe=3 
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U Ye tan^^ cos2(?9e-0 008 2^9^ 


Ml 




f 0.82 0.52 0.22 0.071^ 
y-0.22 0.32 -0.2 0.9 y 


0.85 0.4 1 

-0.33 

0.1 0.22 0.97 


M2 




/ 0.82 0.52 0.22 0.1 \ 
\-0.22 0.1 0.16 0.96/ 


0.85 0.4 1 

0.67 

0.1 0.2 0.96 


M3 




^0.8 0.51 0.22 0.22 ^ 
y0.18 0.25 -0.29 -0.91y 


0.85 0.41 1 

-0.33 

0.1 0.67 0.97 



Table 1: Electron and sterile rows of the three examples of mixing matrices Ml, M2, and 
M3 considered in the text. For each of them we list the value of cos2i}s and the values 
of tan^^ and cos2(?9e — corresponding to two values of the electron fraction Ye- 

Therefore, if ^ 7^ the averaged probability P^^^^"* cannot be equal to P^^^^. Instead, 
from Eq. (I82|) we obtain 

= I cos^ t9/3 [sm\A - sin^ ^ + cos^(?9e - cos^ ^] 

+ sin^ [sin^^^e — cos^ C, + cos^('i9e — sin^ ^] | COS^ Xe cos^ X/3 

N 

+ J2\Uf,k\^\Uek\^. (83) 
k=3 

After simple trigonometric transformations one obtains Eq. fITOl) . 

4 Examples 

In this Section we illustrate the validity of the flavor transition probability fl6T|) and 
the crossing probability flM|) by considering some examples in the simplest case of four- 
neutrino mixing. In this case we have only one sterile neutrino z/^ = Ug^. 

Since we have derived the flavor transition probability without using a specific pa- 
rameterization of the mixing matrix, we continue with this approach. We consider the 



14 



three examples of mixing matrices Ml, M2, and M3 in Tab. [T] In all of them the value 
of de corresponds to the best-fit value of the solar mixing angle obtained in two- and 
three-neutrino mixing fits of solar and KamLAND data [5,78-80]: 

tan^ = tan^ ^sol ^ 0.4 . (84) 

In Tab. [1] we give the electron and sterile rows of the mixing matrix which contain all 
the information needed to calculate the probabilities P^^^y^ and P^^^^,^, from which it is 
immediate to obtain the remaining measurable probability ^^3=^ Pue^Uf^ as 1 — Pu^^u^ — 

In the first example we choose 

Ml: |[/e3|^ = 0.05, |t/e4|^ = 0.005, |f/,ip = 0.05, \Us2? = ^-l- (85) 

We choose relatively large values of |f/sip and |f/s2p in order to have significant z/g — > 
Vs transitions. The values of |f/e3p and |f/e4p are sufficiently small to be compatible 
with the reactor constraints [20,81]. The value of |f/e4p is constrained by the results of 
short-baseline (SBL) reactor experiments in which the effective mixing angle is given by 
sin2 2^9SBL = 4|f/^^|2 _ ^^^^2^ ^^^^ Section 7.7 of Ref. [7]). The Bugey hmit sin^ < 

0.1 for Aml^ > 3 x lO^^eV^ [82] implies that 

|[/e4r< 0.025. (86) 

Neglecting the small value of |?7e4p, the effective mixing angle in long-baseline (LBL) reac- 
tor experiments is given by siia^ 2{}];^^ = 4|?7e3p (1 — l^esP) (see Section 7.8 of Ref. [7]). 
The Chooz limit sin^ 2^9^f ^ < 0.2 for Am|i = Am^^M - 2.5 x lO'^eV^ [83] implies 
that [84,85] 

|f/e3|'< 0.05. (87) 

Therefore, in the first example (Ml) we choose a value of |?7e3p at the level of the 
current experimental upper bound and a value of |?7e4p much smaller that the current 
experimental upper bound. 

The second example has been obtained by choosing 

M2: |[/^3|2 = 0.05, |f/e4p = 0.01, = 0.05 , |[/,2p = 0.01. (88) 

Hence, we consider the same values of |t/e3p and |f/sip as in Ml, but the value of |f/s2p 
is much smaller, and the value of |?7e4p is larger, not far from the current experimental 



upper bound in Eq. fl86|) . 

In the third example we choose 

M3: |f/e3p = 0.05, |[/e4|^ = 0.05, |t/,ip = 0.03, |[/,2|^ = 0.06 . (89) 



The value of |f/e4p is larger than the current experimental upper bound in Eq. fl86l) . It 
is motivated by a possible explanation of the Gallium radioactive source experiments 
anomaly and the MiniBooNE low-energy anomaly through relatively large very-short- 
baseline active-sterile neutrino oscillations [19,21]. This possibility requires a relaxing of 
the reactor constraint on |t/e4p which could be due to an underestimate of systematic 
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Figure 2: Averaged probability of v^. survival as a function of the neutrino energy E 
for the three examples of mixing matrices Ml, M2, and M3 in Tab. [H For each of them 
we consider the two values of the electron fraction Y^. in Tab. [TJ In each plot, the solid 
line is obtained with the analytic expression in Eq. (16T|) . the points are obtained with a 
numerical solution of the evolution equation, the horizontal dashed line shows the value 
of pZ^^„ ■, and the horizontal dotted line shows the value of P^J^^\ 
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Figure 3: Averaged probability of Ue Vg transitions as a function of the neutrino 
energy E for the three examples of mixing matrices Ml, M2, and M3 in Tab. [H For each 
of them we consider the two values of the electron fraction in Tab. [H In each plot, 
the solid line is obtained with the analytic expression in Eq. (1611) . the points are obtained 
with a numerical solution of the evolution equation, the horizontal dashed line shows the 

value of pZ^^u , and the horizontal dotted line shows the value of p']]^!^ . 
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Figure 4: Averaged probability of v^. survival and z/g transitions as functions of the 
neutrino energy E for the three examples of mixing matrices Ml, M2, and M3 in Tab. [T] 
calculated for the BP04 Standard Solar Model density [67]. For each of the three examples 
the lines are obtained with the analytic expression in Eq. fl6T|) and the overlapping points 
are obtained with a numerical solution of the evolution equation. In the left panel we also 
plotted the standard three-neutrino mixing value of Py^^^XE) in the case of negligible 
f/e3 (dash-dotted line). 



uncertainties in the calculation of the reactor neutrino flux or to a violation of the CPT 
symmetry [21]. 

We solved numerically the effective two-neutrino evolution equation fl5^ for the ex- 
ponential electron number density profile 



N^{x) = 245 Na exp (-10.54 x/i?suN) cm" 



(90) 



which is a good approximation of the real electron number density in the solar radiative 
and convective zones (see Ref. [86]). We considered two values of the electron fraction 
Yf.: the realistic value = 0.85, which is a good approximation in the solar radiative 
and convective zones (see Fig. [T]), and the unrealistic value Y^. = 0.1, which is useful to 
test the effects of the neutral-current potential and the case ^ ^ ^e- 

The panels in Figs. [2] and [3] show, respectively, the averaged Ue survival probability and 
the averaged probability of z/g Vg transitions as functions of the neutrino energy E for 
the three examples of mixing matrices Ml, M2, and M3 in Tab. [TJ Notice that in order to 
check beyond any doubt the validity of the analytic flavor transition probability (I^Tj) both 
in the adiabatic and non-adiabatic regimes, we considered a very wide range of energies, 
from 10~^ MeV to 10^ MeV, which is much wider than the energy range of measurable 
solar neutrinos, from about 0.2 MeV to about 15 MeV. Moreover, the inequality fl20l) . 
which allows the approximation in Eq. fl2^ . is not valid for energies not much smaller 
than about 100 MeV, for which the value of Aqq in the center of the Sun is of the same 
order as lAmgj^l ~ Am^rpj^. Hence, we emphasize that the energy values larger than 
those of solar neutrinos in Figs. [2] and [3] are considered only for the purpose of checking 
the validity of the analytic solution of the evolution equation (l36|) . 

From Figs. [2] and [3] one can see that the values of the flavor transition probability 
obtained with the analytic expression (!6T]) and the numerical solution of the evolution 
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equation fl36|l agree perfectly in all the three examples considered. Also, the extreme non- 
adiabatic limit discussed in Section [3] is perfectly described by the analytic approximation. 

The figures illustrate also some different peculiar behavior of the averaged oscillation 
probabilities in the three examples. 

With the realistic electron fraction = 0.85, P^^^^^ has the same qualitative be- 
havior in the three examples, whereas P^^^us is enhanced by the matter effect in Ml 
and M3 and suppressed in M2. This is due to the different signs of cos2'i9s (see Tab. [1]). 
In Ml and M3 cos2'(9s is negative and the adiabatic part of the transition probability 
0.5 cos 2'(9s cos 2'(9g ~ —0.5 cos is positive, leading to an enhancement of P^^^^^. In M2 
the opposite happens. 

In the case of the unrealistic electron fraction Y^. = 0.1, one can see from Figs. [2] and 
[3] that P^J^^l^ 7^ P^^^ve because ^ 7^ '(?e, as discussed in Section |3l The three examples 
have the following interesting differences in the extreme non-adiabatic limit. From Fig. [2] 

one can see that P[,3^e^^ ^ Put-^u^ ^1 ^2 and the inequality is reversed in M3. 
This happens because ^ < in Ml and M2, whereas ^ > in M3. From Fig. [3] one can 

see that P^^^^^^ ^ P^t->vs i^ ^^"^ inequality is reversed in M2 and M3. These 
behaviors are due to the fact that ^ < and cos2'(95 < in Ml, ^ < t^e and cos2't?5 > 
in M2, and ^ > 'de and cos 2^9, < in M3. 

After the check of the analytic expression (1611) for the flavor transition probability in 
the case of the perfect exponential density (1901) and a constant electron number fraction, 
we also checked its validity for the realistic BP04 Standard Solar Model density [67], 
taking into account the variation of the electron number fraction (plotted in Fig. [1]). 
The results of the analytic expression fl6T]) and of the numerical solution of the evolution 
equation fl5B]) are presented in Fig. IHfor the three examples of mixing matrices Ml, M2, 
and M3. We considered only neutrino energies smaller than about 50 MeV, for which 
the inequality (120|) and the approximation (12^ are valid. In this range of energies the 
flavor transition is practically adiabatic and we expect that non-adiabatic effects due to 
the variation of the electron number fraction are negligible. Indeed, from Fig. H] one can 
see that the analytic approximation is very accurate for solar neutrinos and the effects 
of the variation of the electron number fraction, which are neglected in the analytic 
approximation, are really negligible. 

From the left panel in Fig. H] one can see that active-sterile transitions as those in 
the examples Ml, M2, and M3 can change significantly the Ug survival probability with 
respect to the standard three-neutrino mixing value in the case of negligible Ues (dash- 
dotted line). Therefore, it is important to perform an analysis of solar neutrino data [49] 
in the framework of active-sterile neutrino mixing that we have considered in order to 
obtain information on the possible existence of sterile neutrinos and their mixing with 
the active ones. 



5 Conclusions 

In this paper we derived an analytic solution for the flavor transition probabilities of solar 
neutrinos in a general scheme with an arbitrary number of sterile neutrinos, without 
any constraint on the mixing. Hence, we improved the study presented in Ref. [48], 
where it was assumed that the electron neutrino has non-negligible mixing only with the 
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two massive neutrinos which generate the solar squared-mass difference. Here we only 
assumed a realistic hierarchy of neutrino squared-mass differences in which the smallest 
squared-mass difference is effective in solar neutrino oscillations. 

In Section m we have illustrated the validity of the analytical approximation fl^Tl) of the 
flavor transition probability by comparing its value with that obtained with a numerical 
solution of the evolution equation in three examples of the possible mixing matrix in the 
simplest case of four-neutrino mixing. Hence, we are confident that the expression (16T|) 
of the flavor transition probability can be used in a reliable analysis of the solar neutrino 
data which could provide important information on active-sterile neutrino mixing [49]. 
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A Relations Between Bases 

In this paper we used four different bases of the neutrino states and amplitudes: 

1. The flavor basis of amplitudes in Eq. (1121) . related to the flavor states on the left- 
hand side of Eq. ([ID]). 

2. The vacuum mass basis of amplitudes in Eq. (I?li) . related to the vacuum mass basis 
of states in the right-hand side of Eq. (ITUl) . 

3. The effective mass basis in matter of amplitudes in the truncated 1-2 sector, ip'^ 
and ijj^, introduced in Eq. (I57|) . 

4. The effective interaction basis in matter of amplitudes in the truncated 1-2 sector, 
ip\ and ipl, introduced in Eq. (15^ . 

In this Appendix we derive the relations between the effective mass and interaction 
bases in matter and the flavor basis of the amplitudes and the relations of the corre- 
sponding states. 



Since 




TJ* TJ* 
TJ* TJ* 




(91) 



we have 




cos^f/^i + sin^f/^a 
-sin^t/^i + cos^t/^ 



cos^f/;i + sin^f/;2 




(92) 
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Therefore, the effective interaction amphtudes ip\ and ipl for which the matter potential 
in the truncated 1-2 sector is diagonal are linear combinations of the flavor amplitudes 
ipe, ipfi, "ipT, ipsi, ■ ■ ■, i^sMs- Furthermore, from Eqs. fl37|) and flQll) we derive immediately 
the relation between the amplitudes and tp'^ of the effective energy eigenstates in 
matter and the flavor amplitudes ipe,4'ti,'4'T, i^si, ■ ■ ■ , i'sNs '■ 



M 



2 



COS ojUI^ - sin cj[/*2 cos cj[/*^ - sin ujU*^^ 
sin uoU*^ + cos u;f/*2 sin tuf/*^ + cos u;f/*2 



(93) 



In the effective mass and interaction bases in matter the solar neutrino state in Eq. ([7]) 
is given by 



N 



k=l k=3 
2 N 



fc=i 

with the effective massive states in matter 



k=3 



COS LU —smu 
sin uj cos LU 



coscjf/ei — sinci;f/e2 coscuf/^i — smujUi^i2 

sin CUf/el + cos UjUe2 slu CiJ?/^! + COS LjU^2 



and the effective interaction states 



-J) 



cos^ sinA (\v() 
-sin^ cos^y l|z/^) 



cos^M sin^M 
- sin cos 



COS iUe\ + sin ^?7e2 COS + sin 
- sin ^f/ei + cos iUe2 - sin ^f/^i + cos 



Thus, a flavor neutrino states can be expressed as 



N 



k=l 

= (cosuU:, - smujU:2) ) + (smuU:, + cos^f/^^) W^) + J2 U*c.kWk) 

k=3 
N 

= (cosef/^i + sin^U:,) + (- sin^U:, + cos^f/^^) 1^2) + E ^*e.kW^) ' 



TV 



A;=3 



(94) 



(95) 



(96) 



(97) 



for a = e, ^,T, si, . . . , sns- One can see that in general each flavor state contains both 
the effective massive states |z/J^) and |z/^) as well as both the effective interaction states 
|z/{) and 1^1). 
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